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Loop inequalities and confinement 

E.T. Tomboulis'' 

'^Department of Physics, UCLA, Los Angeles, CA 90095-1547, USA 

We consider correlation inequalities that follow from the well-known loop equations of LGT, and their analogues 
in spin systems. They provide a way of bounding long range by short or intermediate range correlations. In several 
cases the method easily reproduces results that otherwise require considerable effort to obtain. In particular, in 
the case of the 2-dimensional 0(N) spin model, where large N analytical results are available, the absence of a 
phase transition and the exponential decay of correlations for all /3 is easily demonstrated. We report on the 
possible application of this technique to the analogous 4-dimensional problem of area law for the Wilson loop in 
LGT at large 



1. Loop equations in LGT 

The loop equations of gauge theory, i.e. the SD 
equations for Wilson loops, are well-known both 
on the lattice and, at least formally, also in the 
continuum m. On the lattice: 






(1) 



Here, and in the following, circles denote gen- 
eral but non-self-intersecting loops. (For self- 
intersecting loops, the l.h.s. contains additional 
terms involving multiple-loop expectations gen- 
erating the infinite sequence of coupled SD loop 
equations. For a simple loop, however, one 
has the closed equation (|l|).) The deformations 
on the r.h.s. involve one plaquette protruding 
indirection^. Also: (3 = l/Ng"^ for U{N), 
f3 = N/[{N^ - l)g^ ] for SU{N). 

2. Basic idea 

In strong coupling expansion for U{N) one 
finds for the 'curly' deformation term 




{2d -3)13- 



(2) 



so the loop equation (|^) gives 



f3 




l-t-/32(2d-3) + 




Now one may iterate any number of times up 
to the number of plaquettes in the minimal loop 
area \A\, and obtain |^: 




< A' 



Max 
all loops 




i.e. area law 




< Const, exp 



-(In-) 1^1 



(4) 



(5) 



provided X 



< 1, 



for 



f3 < 



1 



^ 2{d-l)P 
^ l + /32(2rf-3) 

, . Note that this estimate is com- 
(2d- 3) 

parable to estimates of the convergence radius 
of the strong coupling expansion obtained after 
rather more involved arguments. 

Abstracting from the strong coupling case sug- 
gests the following general approach |2) . Assume 
that the 'curly' term satisfies 




> a(/3) 




(6) 
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for some function a{(3) (for sufficiently large 
loops). Then the loop equation (|^) gives 



(7) 



and area law follows by the above iteration argu- 
ment, previously applied to (||), provided 

W-1) . . 

1 + /3a(/3) < ' 

Attempts to extract area-law from the contin- 
uum loop equations have not been particularly 
successful. Very little seems to have been done 
with the lattice version. Before proceeding to 
our main case of interest, i.e 4-dimensional LGT 
at large /?, we illustrate this general approach in 
some simpler examples. 

3. U(N) LGT in 2 dimensions 

In this case we have the equality: 
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a(/3) 




(9) 



Substituting in the loop equation gives (0) as 
an equality, and, from (^), area law follows if 

a{P) > 2 — — . Now, explicit computation in the 

large N limit gives: 



a{(3) = P, (3<- 
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P 4/3(4/3 - 1) ' ^ 2 



So the above bound on a{P) is satisfied for all 
< /3 < cx). 

4. 0(N) spin model in 2 dimensions 

In the case of the spin model the analog of loop 
equation (|l|) is 

±2 



Ai=±l 



/3 ^ ^ ^ Sn+/1 ■ S,fj 



(10) 
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At=±l 



where S„ denotes an TV-component unit length 
spin at site n. Note that the second term on the 
r.h.s corresponds to the 'curly' term and comes 
with a minus sign. We work in the large N limit. 
To leading order in 1/A^: 



So ([lo| ) becomes 



Sn, ■ St- 



S„-S„)(ll) 



p 



±2 
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+ /3a(/3) E 



(12) 



This is the analog of (mj , and the same reasoning 
implies that ^S„ ■ Smy decays exponentially as 
long as: 

±2 ^ 

a(/3)^ J2 (s„+/,-S„) >4--. (13) 



Now it is a well-known result Q that to leading 



-i- order 



aiP) 



-27r/3 



-I- 



P^OO. (14) 



Hence condition (y_3|) is fulfilled all the way to 
/3 — > oo proving the exponential decay of corre- 
lations and absence of a phase transition in the 
model. Note that this happens by virtue of the 
nonperturbative (exponential) term in (04) sig- 
naling the presence of a spin condensate W. The 
next to leading 1 /N order is rather more involved 
Q with technical complications of the sort one 
encounters in the 4-dim gauge theory case. 

5. 4-dim. SU(N) gauge tiieory at large /3 

We begin by trivially rewriting the curly loop 
in the loop equation (j^) as 

(l-/v)/ \ + K ( \ (15) 





with K an arbitrary parameter to be adjusted 
later for optimizing bounds. 
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Inserting (|l^) in (|l]) , and by a series of manip- 
ulations involving adding and subtracting appro- 
priately chosen loop terms on the r.h.s. and use 
of reflection positivity, one may derive from (|^) 
the inequality: 




P2{d-1) 



- [l + c2G-2cn]i/2 



1 — r(K) 



< 




where c — 1 — k, and r is a somewhat compli- 
cated expression involving differences of ratios of 
loops whose explicit form need not be given here, 
r = 0(1) but is K-dependent. In fact alternative 
versions of this type of inequality with somewhat 
different forms of r may be obtained. Writing for 
the two- and one-plaquette expectations in (p^: 

g = l-E(g), □=1-S](5) 
one has: 

[l + c^[^-2cn] = k2(1-S:) + k2(S-S) 



-(2E - E). 



(17) 



Now expanding for small lattice spacing a and 
using OPE: 

(2E - S) = -a^ ( F^,{x)F'''' [x + a)) + . . . 



^ Ci{a,y)\ + CF.{a,tJL){^g^¥-'{^jL)] 
with leading singular behavior: 
Ci(a, /i) ^ , Cp2 a" . 
We now write the r.h.s. of ( p7| ) as: 

The K independent terms proportional to on 
the r.h.s. of ([T^ ) are uniquely picked out, and 
arise from the gluon condensate: 



(18) 



Oi ~ Const A ^ Const ^ e 



4 „-2bo/s^(M) 



(19) 



After some manipulation the expression multi- 
plying the loop on the l.h.s. in (16) can be written 
as 



(32{d - 1) [ 1 + i 0\'^ + /(g2, /i, «, «) 



where / is given in terms of S, S and the other 
quantities entering in r^n). It involves a series 
in g'^ 'perturbative' part. Note that all depen- 
dence, other than in RG invariant combination 
(as in A), arising through the use of OPE must, 
in principle, cancel among the various terms in /. 

( |l^ ) is optimized by choosing k to maximize 
/. Now since the inequality is rigorously valid, / 
must be negative for all k not equal to the max- 
imizing Kmax, thus rendering the iteration argu- 
ment inapplicable. Otherwise, as it is easily seen, 
unphysical loop behavior would result. For Kmax: 
the optimum that could be achieved is that / van- 
ishes (to within powers of aA) . In such a case ( p^ ) 
would give: 




which may be iterated as above to give area law 
with string tension ~ const A^ . 

Now / at Kmax depends rather delicately on 
the structure of r in (|l6|), and is, unfortunately, 
not easily estimated accurately enough to ascer- 
tain whether it vanishes. This is of course only 
to be expected. The inequality (^6|) must clearly 
be very sharp in order to obtain the right behav- 
ior. As noted there arc different versions of it 
with somewhat different forms of r that may be 
explored, as well as possible ways of successive 
refinement of a given form. These are currently 
under investigation. 
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